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Abstract 

We compute the mid-rapidity densities of pions, kaons, baryons and an- 
tibaryons in Au-Au collisions at y/s = 130 GeV in the Dual Parton Model 
supplemented with final state interaction (comovers interaction). The ra- 
tios B/np ar t (B/np ar t) increase between peripheral (n par t = 18) and central 
(n P art = 350) collisions by a factor 2.4 (2.0) for A's, 4.8 (4.1) for H's and 
16.5 (13.5) for fTs. The ratio K~ /ir~ increases by a factor 1.3 in the same 
centrality range. A comparison with available data is presented. 
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The enhancement of the ratio of yields of strange baryons and antibaryons per 
participant nucleon, observed at CERN-SPS, [|]] @, is one of the main results of 
the Heavy Ion CERN program. A description of these data has been given in Jl in 
the framework of the Dual Parton Model DPM [|J, supplemented with final state 
interaction. The net baryon yield is computed in the same framework taking into 
account the mechanism of baryon stopping, associated with baryon junction transfer 
in rapidity f5|-|9|]. We use its implementation in |§, which describes the SPS data. 

In the absence of nuclear shadowing, the rapidity density of a given type of 



hadron h produced in AA collisions at fixed impact parameter, is given by pfl flT0 



dN AA ~* h 

-(y,b) = n A (b) 



dy 



K,^ (V) + N nJ) (V) + ( 2k - 2 ) N hJ)(v) 

■Is -q a 
h,n(b) 



+ {n{b)-n A {b))2kN^{y) . (1) 



Here 



n{b) = o vv A 2 J d 2 s T A (s) T B (b - s)/a AA (b) (2) 
is the average number of binary collisions and 

n A (b) = A J d 2 s T A (s) [1 - exp(-a pp A T A (b - s)) /cr AA (b) , (3) 

is the average number of participant pairs at fixed impact parameter b. P and T 
denote projectile and target nuclei, k is the average number of inelastic collisions in 
pp and fi(b) = ki>(b) with u(b) = n(b)/n A (b) is the average total number of collisions 



suffered by each nucleon. At ^fs = 130 GeV we have k = 2 JT0| . 

The Nh^(fy(y) in eq. (|1|) are the rapidity distributions of hadron h in each 
individual string. In DPM they are given by convolutions of momentum distribution 
and fragmentation functions^. The first term in (JJ) is the rapidity distribution in 
one NN collision of AA, resulting from the superposition of 2k strings, multiplied 
by the average number of participant pairs. Since in DPM there are two strings 



2 For pions, we use the same fragmentation functions given in p0[ . For simplicity, the same 
form is used for kaons. For pp pair production we take |11| xD"P q (x) = xD^ q (x) ~ (1 — x) 5 and 



xDP(x) — xDP(x) ~ (1 — x) . For the other baryon species an extra a p (0) — a^(0) = 1/2 is added 



to the power of (1 — x) for each strange quark in the baryon |11|. 
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per inelastic collision, the second term, consisting of strings stretched between sea 
quarks and antiquarks, makes up for the total average number of strings 2k n(b). 

It was shown in [HJ that eq. ([!]), supplemented with shadowing corrections, 
leads to values of charged multiplicities at mid-rapidities as a function of centrality 
in agreement with data, both at SPS and RHIC. Here we use the same shadowing 



corrections as in ref. [10] - leading to the lower edge of the shaded area in Fig. 4 of 
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Net Baryon Production . Let us now consider the net baryon production AB = 
B -B. In the standard version of DPM g (or QGSM p[) the leading baryon 
results from the fragmentation of a valence diquark. This component will be called 
diquark preserving (DP). The stopping observed in Pb Pb collisions at SPS has led 
to the introduction of a new mechanism based on the transfer in rapidity of the 
baryon junction |5|-|9|]. Here we follow the formalism in || which describes the SPS 
data. In an AA collision, this component, called diquark-breaking (DB), gives the 
following rapidity distribution of the two net baryons in a single NN collision of AA 



, (y) = c„ ( 6) [zl /2 (i - z + y®-*i* + z l J\i - z_y®-*i* 
d y ) 1/(6) 



(4) 



where Z± = exp(±y — y max ) and u(b) = n(6) /n^(fo). C u ^) is determined from the 
normalization to two at each b. 

The net baryon rapidity distribution in AA collisions is then given by 



dN 



AA^AB 



dy 



-(y,b) = n A {b) 



i/(6) V dy 



(y) 



+ 



i/(6) - 1 



'fJAfAB 



dy 



(y) 



■ (5) 



1/(6) "\"> \ ™» J v(b)_ 

The physical content of eq. (|5p is as follows. Each nucleon interacts in average 
with i/(6) nucleons of the other nucleus. It has been argued in that in only one 
of these collisions the string junction, carrying the baryon number, follows a valence 
diquark, which fragments according to the DP mechanism. In the u(b) — 1 others, 
the string junction is freed from the valence diquark and net baryon production 
takes place according to the DB mechanism, eq. (H). In order to conserve baryon 



number, we have to divide by u(b) and multiply by the number of participating 
nucleons. We obtain in this way eq. (|5|)f]. This equation gives the total net baryon 
density^. 

In order to get the relative densities of each baryon and antibaryon species we 
use simple quark counting rules . Denoting the strangeness suppression factor by 
S/L (with 2L + S — 1), baryons produced out of three sea quarks (which is the case 
for pair production) are given the relative weights 

h = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 . 

for p, n, A + X, H°, S~ and Q, respectively. The various coefficients of ^3 are 
obtained from the power expansion of (2L + S) 3 . In order to take into account the 
decay of E*(1385) into Arc, we redefine the relative rate of A's and S's using the 
empirical rule A = 0.6(S + + £~) - keeping, of course, the total yield of A's plus 
S's unchanged. In this way the normalization constants of all baryon species in 
pair production are determined from one of them. This constant, together with the 
relative normalization of K and 7r, are determined from the data for very peripheral 
collisions. In the calculations we use S = 0.1 (S/L = 0.22). 

For net baryon production two possibilities have been considered. The first one 
is that the behaviour in Z 1 / 2 , eq. (§, is associated to the transfer of the string; 
junction without valence quarks In this case the net baryon is made out of 

three sea quarks and the relevant weights are given by J3. In the second one, eq. 
(f§) is a pre-asymptotic term associated to the transfer of the baryon junction plus 
one valence quark ||. In this case the relevant weights are given by I2, i.e. from 
the various terms in the expansion of (2L + S) 2 . This second possibility is favored 
by the datafj. Since the normalization of the total net baryon yield is determined 
3 In the numerical calculations we neglect the first term of (||) since the DP component gives 
a very small contribution at y* ~ and RHIC energies - about 5% of the DB one for the most 

central bin where its effect is maximal. 

4 In order to conserve strangeness locally, we have added an extra 1/2K + and 1/2K to each 

produced net A (plus S's), an extra K + and K° to each net S and an extra 3/2K + and 3/2K to 

each net fi. 

5 Note, however, that a non-zero value of net omegas has been observed in hA collisions p3j . 
This requires a non- vanishing contribution proportional to I3. However, its effect in AA collisions 
is presumably small since, in this case, the net omegas are almost entirely due to final state 
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from baryon number conservation, there is no extra free normalization constant. 
Moreover, the total net baryon yield is not affected by final state interaction. 



Final State Interactions . The hadronic densities obtained above will be used 
as initial conditions in the gain and loss differential equations which govern final 



state interactions. In the conventional derivation [14] of these equations, one uses 
cylindrical space-time variables and assumes boost invariance. Furthermore, one 
assumes that the dilution in time of the densities is only due to longitudinal motion, 
which leads to a r _1 dependence on the longitudinal proper time r. These equations 
can be written as Jl4| || 

T ~T~ = H a ki Pk <?ik Pi Pk ■ (6) 

dT ki k 

The first term in the r.h.s. of (|6|) describes the production (gain) of particles of 
type i resulting from the interaction of particles k and i. The second term describes 
the loss of particles of type i due to its interaction with particles of type k. In eq. 
© Pi = dNi/dyd 2 s(y,b) are the particles yields per unit rapidity and per unit of 
transverse area, at fixed impact parameter. They can be obtained from the rapidity 
densities ([!]), © using the geometry, i.e. the s-dependence of tia and n, eqs. (fj), (|]). 



The procedure is explained in detail in [[L5| . gm are the corresponding cross-sections 



averaged over the momentum distribution of the colliding particles. 

Equations (|6|) have to be integrated from initial time Tq to freeze-out time Tf. 
They are invariant under the change r — > cr and, thus, the result depends only 
on the ratio T//ro. We use the inverse proportionality between proper time and 
densities and put Ty/r = (dN/dyd 2 s(y,b))/ pf. Here the numerator is given by 
the DPM particles densities. We take Pf — [3/ '-nR^\{dN~ / dy) y *~o = 2 fm~ 2 , which 
corresponds to the charged density per unit rapidity in a pp collision at y/s = 
130 GeV. This density is about 70 % larger [0 than at SPS energies. Since the 
corresponding increase in the AA density is comparable, the average duration time 
of the interaction will be approximately the same at CERN SPS and RHIC - about 
5 to 7 fm. 



interactions. 
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Next, we specify the channels that have been taken into account in our calcula- 
tions. They are 

nN KA(T,) , ttA(S) ^ KE , ttS ^ KQ (7) 
We have also taken into account the strangeness exchange reactions 

ttA(S) <± f^N , ttS ^ KA(E) , vrfi ^ ATS . (8) 

as well as the channels corresponding to ([?]) and (^) for antiparticlesQ We have 
taken = o = 0.2 mb, i.e. a single value for all reactions in ([7]) and (|8]) - the same 
value used in ref. || to describe the CERN SPS data. 



Numerical Results . All our results refer to mid-rapidities. The calculations have 

been performed in the interval —0.35 < y* < 0.35. In Fig. la-Id we show the 

rapidity densities of B, B and B — 5|] versus h~ = dN~ jdr\ = (1 / '1. 17) dN/dy and 

compare them with available data [|T7| — pTQf| . We see that, in first approximation, p, 

p, A and A scale with h~ . Quantitatively, there is a slight decrease with centrality 

of p/h~ and p/h~ ratios, a slight increase of A/h~ and A/h~ and a much larger 

increase for H (E)/h~ and f2 (Q)/h~. In Fig. 2a and 2b we plot the yields of B 

and B per participant normalized to the same ratio for peripheral collisions versus 

n part . The enhancement of B and B increases with the number of strange quarks 

in the baryon. This increase is comparable to the one found at SPS between pA 

and central PbPb collisions - somewhat larger for antibaryons. The ratio K~/tt~ 

increases by 30 % in the same centrality range, between 1.1 and 1.4 in agreement 

with present data [|17]. The ratios B/B have a mild decrease with centrality of 

6 To be precise, of all possible charge combinations in reactions (0), we have only kept those 
involving the annihilation of a light q-q pair and production of an s-s in the s-channel. The other 
reactions, involving three quarks in the i-channel intermediate state, have substantially smaller 
cross-sections and have been neglected. All channels involving 7r° have been taken with cross- 
section a/2 since only one of the uu and dd components of ir° can participate to a given charge 

combination. For details see the second paper of ||. 

7 A Monte Carlo calculation in a similar framework with string fusion can be found in |j~6|] . A 

net proton rapidity density of about 10 for central Au Au collisions at mid-rapidities at RHIC was 

first predicted in S using a stopping mechanism similar to the one considered here. 
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about 15 % for all baryon species - which is also seen in the data |2T|. Our values 
for N ch /N± x = 1/2 are : p/p = 0.69, A/A = 0.72, E/S = 0.79, Q/H = 0.83Q, to be 
compared with the measured values [[Tj| : 

p/p = 0.63 ±0.02 ±0.06 , A/A = 0.73 ± 0.03 , H/S = 0.83 ± 0.03 ± 0.05 . 

The ratio K + /K~ = 1.1 and has a mild increase with centrality, a feature also seen 
in the data. 

As explained above only one parameter has been adjusted in order to determine 



the absolute yields of baryons and antibaryons. Since both the PHENIX [17| and 
the STAR data are not corrected for feed-down from weak decays, this free 
parameter has to be re-determined after these corrections are known. It will then 
be possible to compare its value with the one obtained from other sets of data, in 
particular pp. Likewise, the corrected yields of net protons will determine the exact 
amount of stopping and should allow to decide whether or not there is "anomalous" 
stopping in AA, i.e. an excess as compared to an extrapolation from pp and pA. At 
present there is no clear indication of anomalous stopping. 



After completion of this work, the PHENIX collaboration |24[] has published the 
yields of p and p for the 5 % most central events corrected for feed-down. The cor- 
rections are of 30 %.The systematic error is 20 %.Also, C.Roy (STAR collaboration) 
communicated to us a recent STAR preliminary result on the mid-rapidity densities 
of S~ and H + for the 14 % most central events. These values are 60 to 70 % higher 
than our results. Their systematic error is 20 %. 

Physical interpretation . Before final state interactions, all ratios K/h~, B/h~ 
and B/h~ decrease slightly with increasing centrality. This effect is rather marginal 
at RHIC energies and mid-rapidities. 

The final state interactions (0), @ lead to a gain of strange particle yields. The 
reason for this is the following. In the first direct reaction (0) we have p n > px, 
Pn > Pa, P-kPn PkPa- The same is true for all direct reaction (|7|). In view of 
that, the effect of the inverse reactions (0) is small. On the contrary, in all reactions 
8 In ref. p], the relative weights of net baryons were given by the factors 0.5(/2 -I-/3) - instead of 
J 2 - In this case the values of the ratios are p/p = 0.70, A/A = 0.71, 3/3 = 0.76 and 0/f2 = 0.78. 
Their increase with the number of strange quarks in the baryon is smaller. 
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([8]), the product of densities in the initial and final state are comparable and the 
direct and inverse reactions tend to compensate with each other. Baryons with 
the largest strange quark content, which find themselves at the end of the chain of 
direct reactions (0) and have the smallest yield before final state interaction, have 
the largest enhancement. Moreover, the gain in the yield of strange baryons is larger 
than the one of antibaryons since Pb > Pb- Furthermore, the enhancement of all 
baryon species increases with centrality, since the gain, resulting from the first term 
in eq. (|6|), contains a product of densities and thus, increases quadratically with 
increasing centrality. 

Although the inverse slopes ("temperature") have not been discussed here, let us 
note that in DPM they are approximately the same for all baryons and antibaryons 
both before and after final state interaction - the effect of final state interaction 
being rather small [|16j [ [23f . 

In conclusion, we would like to emphasize the fact that in DPM (before final 
state interaction) the rapidity density of charged particle per participant increases 
with centrality. This increase is larger for low centralities ||10|| . This has an im- 
portant effect on both the size and the pattern of strangeness enhancement in our 
results. It explains why the departure from a linear increase of S's and JTs (concave 
shape) seen in Figs, lc and Id is also more pronounced for lower centralities. It 
leads to the convexity in the centrality dependence of the yields of hyperons and an- 
tihyperons per participant in Figs. 2 (Note, however, a change of curvature for very 
peripheral collisions where the effect of final state interaction is negligible). This 
centrality pattern is a distinctive feature as well as a firm prediction of our approach. 
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Figure Caption : 

Figure 1. (a) Calculated values of rapidity densities of p (solid line), p (dashed 
line), and p — p (dotted line) at mid rapidities, \y*\ < 0.35, are plotted as a function 



of dN I dr], and compared with PHENIX data [17] ; (b) same for A and A compared 



to preliminary STAR data |HJ ; (c) same for E and H ; (d) same for Q and f2. 

Figure 2. Calculated values of the ratios B/n part (a) and B/n part (b), normalized 
to the same ratio for peripheral collisions (n part = 18), are plotted as a function of 
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Fig. 1 
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Fig. 2.b 
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The enhancement of the ratio of yields of strange baryons and antibaryons per 
participant nucleon, observed at CERN-SPS, [1] [2], is one of the main results of 
the Heavy Ion CERN program. A description of these data has been given in [3] in 
the framework of the Dual Parton Model DPM [4], supplemented with final state 
interaction. The net baryon yield is computed in the same framework taking into 
account the mechanism of baryon stopping, associated with baryon junction transfer 
in rapidity [5-9]. We use its implementation in [3], which describes the SPS data. 

In the absence of nuclear shadowing, the rapidity density of a given type of 
hadron h produced in AA collisions at fixed impact parameter, is given by [4] [10] 



dN AA ~" h 

-(y,b) = n A (b) 



dy 



<(/(?/) + <( 6 f (y) + (2* " 2) JVQjM 
+ (n(b)-n A (b))2kN*-fi(y) . (1) 



Here 



n(b) = a pp A 2 J d 2 s T A (s) T B (b - s)/a AA {b) (2) 
is the average number of binary collisions and 

n A (b) = A J d 2 s T A (s) [1 - exp(-a pp A T A (b - s)} /a AA {b) , (3) 

is the average number of participant pairs at fixed impact parameter b. P and T 

denote projectile and target nuclei, k is the average number of inelastic collisions in 

pp and ji(b) = kv(b) with v(b) = n(b)/n A (b) is the average total number of collisions 

suffered by each nucleon. At y/s = 130 GeV we have k = 2 [10]. 

The N h ^ b )(y) in eq. (1) are the rapidity distributions of hadron h in each 

individual string. In DPM they are given by convolutions of momentum distribution 

and fragmentation functions 2 . The first term in (1) is the rapidity distribution in 

one NN collision of AA, resulting from the superposition of 2k strings, multiplied 

by the average number of participant pairs. Since in DPM there are two strings 

2 For pions, we use the same fragmentation functions given in [10]. For simplicity, the same 
form is used for kaons. For pp pair production we take [11] xD^ q (x) = xDP q (x) ~ (1 — a;) 5 and 
xDv(x) = xDv(x) ~ (1 — x) 3 . For the other baryon species an extra 0^(0) — a^(0) = 1/2 is added 
to the power of (1 — x) for each strange quark in the baryon [11]. 
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per inelastic collision, the second term, consisting of strings stretched between sea 
quarks and antiquaries, makes up for the total average number of strings 2k n(b). 

It was shown in [10] that eq. (1), supplemented with shadowing corrections, 
leads to values of charged multiplicities at mid-rapidities as a function of centrality 
in agreement with data, both at SPS and RHIC. Here we use the same shadowing 
corrections as in ref. [10] - leading to the lower edge of the shaded area in Fig. 4 of 
[10]. 

Net Baryon Production . Let us now consider the net baryon production AB = 
B -B. In the standard version of DPM [4] (or QGSM [12]) the leading baryon 
results from the fragmentation of a valence diquark. This component will be called 
diquark preserving (DP). The stopping observed in Pb Pb collisions at SPS has led 
to the introduction of a new mechanism based on the transfer in rapidity of the 
baryon junction [5-9]. Here we follow the formalism in [3] which describes the SPS 
data. In an AA collision, this component, called diquark-breaking (DB), gives the 
following rapidity distribution of the two net baryons in a single NN collision of AA 
[3] 



dy 



(y)) = c v[b) \z l '\i - z + y(w + z l J 2 (i - z_)^)- 3 / 2 ] (4) 

/ 1/(6) 



where Z± = exp(±y — y max ) and v(b) = n(b)/riA{b). C v (p) is determined from the 
normalization to two at each b. 

The net baryon rapidity distribution in AA collisions is then given by 



dN AA^AB 

dy 



(y,b) = n A (b) 



-(b) { dy [y) )„ {b) + i/(6) 1. dy [y) )„ (b) 



(5) 



The physical content of eq. (5) is as follows. Each nucleon interacts in average 
with v{b) nucleons of the other nucleus. It has been argued in [3] that in only one 
of these collisions the string junction, carrying the baryon number, follows a valence 
diquark, which fragments according to the DP mechanism. In the v(b) — 1 others, 
the string junction is freed from the valence diquark and net baryon production 
takes place according to the DB mechanism, eq. (4). In order to conserve baryon 



number, we have to divide by v(b) and multiply by the number of participating 
nucleons. We obtain in this way eq. (5) 3 . This equation gives the total net baryon 
density 4 . 

In order to get the relative densities of each baryon and antibaryon species we 
use simple quark counting rules [3]. Denoting the strangeness suppression factor by 
S/L (with 2L + S = 1), baryons produced out of three sea quarks (which is the case 
for pair production) are given the relative weights 

I 3 = 4L 3 : 4L 3 : 12L 2 S : 3LS 2 : 3LS 2 : S 3 . 

for p, n, A + X, S°, S~ and Q, respectively. The various coefficients of I 3 are 
obtained from the power expansion of (2L + S) 3 . In order to take into account the 
decay of E*(1385) into Air, we redefine the relative rate of A's and E's using the 
empirical rule A = 0.6(S + + E~) - keeping, of course, the total yield of A's plus 
E's unchanged. In this way the normalization constants of all baryon species in 
pair production are determined from one of them. This constant, together with the 
relative normalization of K and n, are determined from the data for very peripheral 
collisions. In the calculations we use S = 0.1 (S/L = 0.22). 

For net baryon production two possibilities have been considered. The first one 
is that the behaviour in Z 1 / 2 , eq. (4), is associated to the transfer of the string 
junction without valence quarks [5,7,9]. In this case the net baryon is made out of 
three sea quarks and the relevant weights are given by I 3 . In the second one, eq. 
(4) is a pre-asymptotic term associated to the transfer of the baryon junction plus 
one valence quark [6]. In this case the relevant weights are given by I2, i.e. from 
the various terms in the expansion of (2L + S) 2 . This second possibility is favored 
by the data 5 . Since the normalization of the total net baryon yield is determined 
3 In the numerical calculations we neglect the first term of (5) since the DP component gives 
a very small contribution at y* ~ and RHIC energies - about 5% of the DB one for the most 

central bin where its effect is maximal. 

4 In order to conserve strangeness locally, we have added an extra 1/2K + and 1/2K to each 

produced net A (plus S's), an extra K + and K° to each net H and an extra 3/2K + and 3/2K to 

each net ft. 

5 Note, however, that a non-zero value of net omegas has been observed in hA collisions [13]. 
This requires a non- vanishing contribution proportional to ^3. However, its effect in AA collisions 
is presumably small since, in this case, the net omegas are almost entirely due to final state 
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from baryon number conservation, there is no extra free normalization constant. 
Moreover, the total net baryon yield is not affected by final state interaction. 

Final State Interactions . The hadronic densities obtained above will be used 
as initial conditions in the gain and loss differential equations which govern final 
state interactions. In the conventional derivation [14] of these equations, one uses 
cylindrical space-time variables and assumes boost invariance. Furthermore, one 
assumes that the dilution in time of the densities is only due to longitudinal motion, 
which leads to a r _1 dependence on the longitudinal proper time r. These equations 
can be written as [14] [3] 

r^ 1 = ^2 ou Pk Pi~Yj a ik Pi Pk ■ (6) 

ClT kl k 

The first term in the r.h.s. of (6) describes the production (gain) of particles of 
type i resulting from the interaction of particles k and I. The second term describes 
the loss of particles of type i due to its interaction with particles of type k. In eq. 
(6) pi = dNi/dyd 2 s(y,b) are the particles yields per unit rapidity and per unit of 
transverse area, at fixed impact parameter. They can be obtained from the rapidity 
densities (1), (5) using the geometry, i.e. the s-dependence of ha and n, eqs. (2), (3). 
The procedure is explained in detail in [15]. cry are the corresponding cross-sections 
averaged over the momentum distribution of the colliding particles. 

Equations (6) have to be integrated from initial time To to freeze-out time tj. 
They are invariant under the change r — >■ cr and, thus, the result depends only 
on the ratio tj/tq. We use the inverse proportionality between proper time and 
densities and put tj/tq = (dN/dyd 2 s(y,b))/ pf. Here the numerator is given by 
the DPM particles densities. We take pf = [3/7rR 2 ](dN~ /dy) y *^ = 2 fm~ 2 , which 
corresponds to the charged density per unit rapidity in a pp collision at y/s = 
130 GeV. This density is about 70 % larger [10] than at SPS energies. Since the 
corresponding increase in the AA density is comparable, the average duration time 
of the interaction will be approximately the same at CERN SPS and RHIC - about 
5 to 7 fm. 
interactions. 
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Next, we specify the channels that have been taken into account in our calcula- 
tions. They are 



irN KA(S) , ttA(£) KE 



(7) 



We have also taken into account the strangeness exchange reactions 



ttA(E) KN , ttS ^ KA(E) , irQ^KZ . 



(8) 



as well as the channels corresponding to (7) and (8) for antiparticles 6 . We have 
taken = a = 0.2 mb, i.e. a single value for all reactions in (7) and (8) - the same 
value used in ref. [3] to describe the CERN SPS data. 

Numerical Results . All our results refer to mid-rapidities. The calculations have 

been performed in the interval —0.35 < y* < 0.35. In Fig. la-Id we show the 

rapidity densities of B, B and B — B 7 versus h~ = dN~ jdr\ = (1 / 1. 17) dN/dy and 

compare them with available data [17-19]. We see that, in first approximation, p, 

p, A and A scale with h~ . Quantitatively, there is a slight decrease with centrality 

of pjhr and pjhr ratios, a slight increase of kjhr and kjhr and a much larger 

increase for H and Q (Q)/h~. In Fig. 2a and 2b we plot the yields of B 

and B per participant normalized to the same ratio for peripheral collisions versus 

n part- The enhancement of B and B increases with the number of strange quarks 

in the baryon. This increase is comparable to the one found at SPS between pA 

and central PbPb collisions - somewhat larger for antibaryons. The ratio K~ /tt~ 

increases by 30 % in the same centrality range, between 1.1 and 1.4 in agreement 

with present data [17]. The ratios B/B have a mild decrease with centrality of 

6 To be precise, of all possible charge combinations in reactions (7), we have only kept those 
involving the annihilation of a light q-q pair and production of an s-s in the s-channel. The other 
reactions, involving three quarks in the i-channel intermediate state, have substantially smaller 
cross-sections and have been neglected. All channels involving ir° have been taken with cross- 
section cr/2 since only one of the uu and dd components of 7r° can participate to a given charge 

combination. For details see the second paper of [3]. 

7 A Monte Carlo calculation in a similar framework with string fusion can be found in [16]. A 

net proton rapidity density of about 10 for central Au Au collisions at mid-rapidities at RHIC was 

first predicted in [8] using a stopping mechanism similar to the one considered here. 
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about 15 % for all baryon species - which is also seen in the data [21]. Our values 
for N ch /N^ ax = 1/2 are : p/p = 0.69, A/A = 0.72, H/H = 0.79, tl/U = 0.83 8 , to be 
compared with the measured values [18] : 

p/p = 0.63 ±0.02 ±0.06 , A/A = 0.73 ± 0.03 , S/S = 0.83 ± 0.03 ± 0.05 . 

The ratio K + / K~ = 1.1 and has a mild increase with centrality, a feature also seen 
in the data. 

As explained above only one parameter has been adjusted in order to determine 
the absolute yields of baryons and antibaryons. Since both the PHENIX [17] and 
the STAR data [19] are not corrected for feed-down from weak decays, this free 
parameter has to be re-determined after these corrections are known. It will then 
be possible to compare its value with the one obtained from other sets of data, in 
particular pp. Likewise, the corrected yields of net protons will determine the exact 
amount of stopping and should allow to decide whether or not there is "anomalous" 
stopping in AA, i.e. an excess as compared to an extrapolation from pp and pA. At 
present there is no clear indication of anomalous stopping. 

After completion of this work, the PHENIX collaboration [24] has published the 
yields of p and p for the 5 % most central events corrected for feed-down. The cor- 
rections are of 30 %.The systematic error is 20 %.Also, C.Roy (STAR collaboration) 
communicated to us a recent STAR preliminary result on the mid-rapidity densities 
of £~ and S + for the 14 % most central events. These values are 60 to 70 % higher 
than our results. Their systematic error is 20 %. 

Physical interpretation . Before final state interactions, all ratios K/h~, B/h~ 
and B/h~ decrease slightly with increasing centrality. This effect is rather marginal 
at RHIC energies and mid-rapidities. 

The final state interactions (7), (8) lead to a gain of strange particle yields. The 
reason for this is the following. In the first direct reaction (7) we have p v > px, 
Pn > Pa, P-kPn 3> PkPk- The same is true for all direct reaction (7). In view of 
that, the effect of the inverse reactions (7) is small. On the contrary, in all reactions 
8 In ref. [3], the relative weights of net baryons were given by the factors 0.5(/2 +^3) - instead of 
h- In this case the values of the ratios are p/p = 0.70, A/A = 0.71, S/S = 0.76 and O/O = 0.78. 
Their increase with the number of strange quarks in the baryon is smaller. 
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(8), the product of densities in the initial and final state are comparable and the 
direct and inverse reactions tend to compensate with each other. Baryons with 
the largest strange quark content, which find themselves at the end of the chain of 
direct reactions (7) and have the smallest yield before final state interaction, have 
the largest enhancement. Moreover, the gain in the yield of strange baryons is larger 
than the one of antibaryons since p# > p^. Furthermore, the enhancement of all 
baryon species increases with centrality, since the gain, resulting from the first term 
in eq. (6), contains a product of densities and thus, increases quadratically with 
increasing centrality. 

Although the inverse slopes ("temperature") have not been discussed here, let us 
note that in DPM they are approximately the same for all baryons and antibaryons 
both before and after final state interaction - the effect of final state interaction 
being rather small [16] [23]. 

In conclusion, we would like to emphasize the fact that in DPM (before final 
state interaction) the rapidity density of charged particle per participant increases 
with centrality. This increase is larger for low centralities [10]. This has an im- 
portant effect on both the size and the pattern of strangeness enhancement in our 
results. It explains why the departure from a linear increase of S's and fi's (concave 
shape) seen in Figs, lc and Id is also more pronounced for lower centralities. It 
leads to the convexity in the centrality dependence of the yields of hyperons and an- 
tihyperons per participant in Figs. 2 (Note, however, a change of curvature for very 
peripheral collisions where the effect of final state interaction is negligible). This 
centrality pattern is a distinctive feature as well as a firm prediction of our approach. 
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Figure Caption : 

Figure 1. (a) Calculated values of rapidity densities of p (solid line), p (dashed 
line), and p — p (dotted line) at mid rapidities, \y*\ < 0.35, are plotted as a function 
of dN~ /dr], and compared with PHENIX data [17] ; (b) same for A and A compared 
to preliminary STAR data [19] ; (c) same for S~ and S + ; (d) same for and fi. 

Figure 2. Calculated values of the ratios B/n part (a) and B/n part (b), normalized 
to the same ratio for peripheral collisions (n part = 18), are plotted as a function of 

Impart ■ 
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Fig. 1 




Fig. 2.a 



Q 



A 

P 



50 100 150 200 250 300 



350 400 
n part 



12 



Fig. 2.b 
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